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SUMMARY 


The formation of soot always accompanies the combustion of hydrocarbon 
fuels. Sooting causes several undesirable effects and becomes particularly 
severe with synthetic fuels. Although the recent "oil crisis" seems to be over, 
the development and future use of synthetic fuels are inevitable. Soot 
formation is one of the challenging problems in the development of these new 
fuels . 

The Combustion Laboratory at Louisiana State University has long been 
engaged in studies of soot formation. Both experimental and theoretical studies 
have been conducted. The ultimate goal of this program is__the~ development of a 
detailed kinetic model for soot formation during the combustion of hydrocarbons. 
The objective of the experimental part of the contract is to establish 
quantitative relationships between the parameters of soot formation (i.e. soot 
yields, induction times, and rates of formation) and experimental conditions, 
such as temperature, pressure, composition of the mixture and molecular 
structure of the fuel molecules. The experiments have involved pyrolysis and 
oxidation of a variety of hydrocarbon fuels using state-of-the-art shock-tube, 
laser diagnostic and data processing techniques. 

As a result of our activity in recent years, a conceptual model for soot 
formation during the pyrolysis of aromatic hydrocarbons was developed. Under 
the current contract, a conceptual model for soot formation during pyrolysis of 
non-aromatic hydrocarbons was completed. It was suggested that the incipient 
soot formation from hydrocarbons must follow the route of consecutive production 
of conjugated reactive structures. The difference in soot formation 
characteristics between various hydrocarbons is determined by the initiation 
process, i.e. by the reactions leading to these reactive structures. 
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Soot formation in toluene-, benzene-, and acetylene-oxygen-argon mixtures 
was investigated to study soot formation in a combustion environment. It was 
observed that high concentrations of oxygen completely suppress soot formation. 
The addition of oxygen at relatively low concentrations uniformly suppresses 
soot formation at high pressures, while at relatively lower pressures it 
suppresses soot formation at higher temperatures while promoting soot production 
at lower temperatures. The observed behavior may indicate that oxidation 
reactions compete with ring fragmentation. The main conclusion to be drawn from 
the results of this work is that the soot formation mechanism is probably the 
same for the pyrolysis and oxidation of hydrocarbons. That is, the addition of 
oxygen does not alter the soot route but rather promotes or inhibits this route 
by means of competitive reactions. 



3 


INTRODUCTION 


The formation of soot always accompanies the combustion of hydrocarbon 
fuels. Soot production causes several undesirable effects such as air 
pollution, emission of carcinogenic materials and a decrease in the efficiency 
of practical combustion systems. Sooting is particularly severe with synthetic 
fuels due to their higher aromatic content. 

The subject of soot formation attracted particular attention with the 
recent "oil crisis". Although the "crisis" seems to be over, the development 
and future use of synthetic fuels are inevitable. One must realize that the 
many scientific and technological problems associated with future fuels may 
require time to be resolved. Soot formation is certainly one of those 
challenging problems. 

The Combustion Laboratory at Louisiana State University has long been 
engaged in studies of soot formation. Both experimental and theoretical studies 
have been conducted. The experiments have involved pyrolysis and oxidation of a 
variety of fuels using state-of-the art shock tube, laser diagnostic and data 
processing techniques. The ultimate goal of this program is the development of 
a detailed kinetic model for soot formation during the combustion of 
hydrocarbons . 
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1 . Literature Survey 

During the last few years a number of excellent reviews have appeared 
covering various aspects of soot formation [ 1 — 1 0 ] . Here we will concentrate on 
the subject of the ongoing research program, namely shock-tube oxidation of 
hydrocarbons and closely related topics. 

a) Oxidation and Pyrolysis of Aromatic Hydrocarbons in Shock Tubes . 

Early shock tube works were not concerned with formation of soot; rather, 
they concentrated on ignition delay times and the initial stages of pyrolysis 
and oxidation of aromatic fuels. Kogarko and Borisov [ll] measured the 
induction periods for 3% benzene - air mixtures. Miyama extended their 
measurements to a range of benzene concentrations [12]. He also determined 
ignition delays for other aromatic hydrocarbons: toluene, o-xylene, m-xylene, 
p-xylene, ethylbenzene, propylbenzene and 1 ,3,5-trimethylbenzene [13]. The most 
recent measurements of ignition delay times for benzene and toluene were 
reported by Matula and Farmer [14], Burcat et al. [15], and Durgaprasad [16]. 
Table 1 summarizes the ignition delay data, where I is the ignition delay time 
defined as the time from the onset of reaction to the instant of ignition. The 
bracketed quantities in Table I, [fuel], [oxygen], [inert] are molar 
concentrations (mol/cm 3 ). 



TABLE I. IGNITION DELAY TIMES 
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The first shock tube studies with aromatic hydrocarbons were probably done 
in 1963. Bauer and Aten [17] measured absorption spectra of benzene and of 
hexaf luorobenzene diluted in argon over the temperature range from 690 to 1900 
K. They suggested a chain mechanism for the initial stages of benzene 
pyrolysis which can be abbreviated to 


C 6 H 6 


P 2 + P 3 + P 4 


where: 


C 4 H 2 


C 12 H 10 

C 8 H 6 


C 10 H 8 ' 


Orr [18] investigated combustion of dilute fuel-oxygen-argon mixtures for 
ethylene, acetylene, isooctane, n-heptane, benzene, and ^-CO. The rates of 
formation of reactants, intermediates and products were studied by observing the 
infrared emissions. Orr found that production of CO from benzene was 
considerably slower than from other fuels studied. 

Asaba and Fujii [19] investigated pyrolysis of 10-20% benzene-argon 
mixtures at temperatures from 1400 to 1900 K. A single-pulse shock-tube 
technique accompanied by light absorption was employed. They found that the 
pyrolysis reaction was slowed considerably by the addition of hydrogen and 
promoted by the addition of methane. To explain the experimental results, the 
authors proposed a new chain mechanism for benzene pyrolysis: 



C 2 H 2 + C 4 H 3 


-*■ Soot 


Later, Fujii and Asaba [20] and Fujii et al. [21] extended this study to include 
the oxidation of rich and lean mixtures of benzene and oxygen. Their reaction 
mechanism provided a reasonable description of the ignition characteristics 
measured by Miyama [13] but it failed to predict the correct concentrations for 
all the species. McLain et al. [22] modified the mechanism of Fujii and Asaba 
for benzene oxidation and suggested a reaction mechanism for toluene oxidation. 
The latter mechanisms predicted the ignition delays reasonably well but over- 
estimated the concentration of CO and underestimated the concentration of CO^. 

Troe and co-workers [23] reported rate constants for the unimolecular 
decomposition of toluene, C^H^CH^, and of the benzyl radical, CgH^-CI^. They 
clearly demonstrated that both decompositions are within the fall-off region at 
pressures up to 30 atm. Dzarnoski et al. [24] measured the rate constant of the 
reaction 

CH 3 +C 6 H 5 CH 3 * CH 4+ C 6 H 5 CH 2 . 

A number of rate constants for reactions of aromatic hydrocarbons were recently 
determined by methods other than shock tube [25-31]. 

Kern et al. [32] have recently investigated the product distribution in the 
shock tube pyrolysis of toluene and benzene. The major products detected by a 
time-of-flight mass spectrometer were acetylene and polyacetylenes. No products 
with a mass greater than 92 (for toluene) or 78 (for benzene) were detected. 
Their results parallel those of Smith [33,34], which were obtained in high- 
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temperature-low-pressure-Knudsen-cell experiments. The results of very-low- 
pressure pyrolyses of several aromatic species have recently appeared in the 
literature [35-37]. A new mechanism for high-temperature oxidation of 
aromatic hydrocarbons has been recently suggested by Glassman and co-workers 
[38]. 

Mar'yasin and Nabutovskii [39,40] were probably the first to include 
quantitative, although very approximate, measurements of soot in the single- 
pulse shock-tube studies of the pyrolysis of aromatic hydrocarbons. They 
pyrolyzed 3-5% benzene-argon mixtures behind reflected shock waves at 
temperatures from 1400 to 2500 K. Concentrations of hydrogen, methane, 
ethylene, acetylene, vinylacetylene, diacetylene and benzene were analyzed by 
gas-chromatograph and the amount of soot formed was determined by weighing the 
soot deposited on the walls and optical windows of the shock tube. Induction 
times for soot appearance were quantified and correlated to the initial 
temperature and benzene concentration. 

Park and Appleton [41 ] employed a laser extinction technique (633 nm, 
He-Ne) to study soot oxidation rates in an oxygen atmosphere behind reflected 
shock waves at temperatures of 1700-4000 K and pressures of 0.05 - 13 atm. 

In their classical work, Graham et al. [42,43] investigated soot formation 
in pyrolyses of benzene, toluene, ethylbenzene and other aromatic fuels highly 
diluted in argon behind incident shock waves. Their experimental conditions 
covered a temperature range of 1600-2300 K at a total carbon concentration in 
the shock-heated gases of (1. 5-3.0) x 10 atoms/cm . Soot was monitored by 
laser-light extinction measurements at three different wavelengths (488 nm of 
argon- ion laser and 632.8 nm and 3.39 pm of He-Ne laser) and by light scattering 
(488 nm of argon-ion laser) . Samples of soot particles formed in the incident 
shock flows were examined by electron microscope. Their main conclusions were: 
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a) Soot yield passes through a well-pronounced maximum around 1750 

K. Complete conversion of fuel carbon to soot was assumed at the 
maximum point. 

b) "... aromatic hydrocarbons can form soot by two very different 
pathways, one direct and one indirect." That is: 

fast soot 

(via direct 

route) 

soot 

CLft slow (via indirect 

C!TH X route) 

3 x 

Thus, the maximum in soot yield was explained by competition 
between condensation and fragmentation of the aromatic ring. 

c) The authors suggested a model of "f ree-molecule coagulation of an 
aerosol comprised of spherical sticky particles whose size 
distribution is of the self-preserving form". Interpreting the 
experimental results in light of this model, they estimated 

that the particle did not grow larger then about 40 nm. The 
later development of this aspect is reported in [44,45]. 

Wang et al. [46] conducted experiments similar to those of Graham et al. 
[42] but behind reflected shock waves. They confirmed the existence of a 
maximum for soot production as a function of temperature. The maximum 
conversion of carbon atoms to soot was estimated to be approximately 80%. Wang 
et al. [46,47] also determined: 

a) hydrogen and oxygen suppress soot formation; 

b) pressure does not have significant effect on soot formation; 

c) the amount of soot formed increases with increasing initial 
concentration of aromatic hydrocarbon; 
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d) induction times of soot appearance were determined as a function 
of the temperature and initial composition of the mixture; 

e) based on thermodynamic considerations, the precursor to soot was 
suggested to be C 96 H 24- 

Vaughn et al. [48,49] conducted a single-pulse shock- tube study of thermal 
decomposition of benzene diluted in argon. The reaction was studied over a 
temperature range of 1300-2300 K with initial benzene concentrations of (0.4- 
1.25) x 10 6 mole/cm and reaction times of 0.1 to 3.0 msec. The stable products 
formed were analyzed by gas chromatography and mass spectrometry. The amount of 
soot formed was determined by gravimetric methods . Product analysis showed that 
the principal products were acetylene and styrene. Smaller concentrations of 
diacetylene, methane, vinylacetylene , ethylene and toluene were observed. Their 
main conclusions were: 

a) acetylene is principally formed directly from ring scission, 
compounds are formed primarily via acetylene recombination, and 
styrene is formed mainly by an "acetylenic reaction with the benzene 
ring" ; 

b) soot yields increase from zero at 1300 K to a plateau value of 0.8 at 
about 1900 K; 

c) the amount of soot formed decreases when the initial concentration of 
aromatic hydrocarbon is increased. (Note that this observation is 
opposite to that of Wang et al. [46]). 

Later, Nelson et al. [50] extended these studies to include laser extinction 
measurements. Soot yields determined by the optical method exhibited 
maxima similar to those of Graham et al. [42] and Wang et al. [46]. However, 
the dependency on the initial concentration of benzene obtained was similar to 



12 


that of Vaughn et al. [49]. Nelson et al. [50] reported that soot yields formed 
in pyrolyses of toluene were larger than those obtained in pyrolyses of benzene. 
In their work induction times for soot appearance were determined by a laser 
extinction technique and soot samples were analyzed by a transmission electron 
microscope. The authors also suggested that "the primary impact of ring 
fragmentation is not on the nucleation of soot, but rather on the surface 
growth" . 

Evans and Williams [51] characterized soot particles formed behind 
reflected shocks by a transmission electron microscope. They found that the 
agglomeration of soot particles depends on experimental conditions. Williams' 
graduate student, Durgaprasad [16], reported that the pyrolysis of toluene 
resulted in larger amounts of soot than that of benzene. Farmer, Edelman, and 
co-workers [52,53] attempted quasi-global modeling to predict soot emissions in 
practical combustion environments. 

Frenklach et al. [54] re-investigated soot formation in the pyrolysis of 
toluene behind reflected shock waves at conditions similar to those of Wang et 
al. [46,47] and Kern et al. [32]. Soot was monitored by attenuation of a laser 
beam in the visible (632.8 nm) and the infrared (3.39 pm) regions of the 
spectrum. The experimental data indicated that the position of measured maximum 
in soot yield is not universal, but rather is dependent on controllable 
experimental variables including observation time, total pressure and initial 
hydrocarbon concentration, as well as the wavelength employed in the 
measurements. The observed pressure effect could not be rationalized within 
Graham's model which was discussed earlier in the text. Frenklach et al. [54] 
proposed a new conceptual model for soot formation during the high temperature 
pyrolysis of aromatic hydrocarbons. The model is described in terms of the 


kinetic skeleton: 
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A 

A + X 





X 

s 


(Rl) 

(R2) 


where A, X, and S denote an intact aromatic ring, intermediate species and soot, 
respectively. Reactions (Rl) and (R2) generalize the fragmentation of the ring 
and the free-radical polymerization, respectively. This new model explains and 
unifies the various experimental facts and theoretical hypotheses. Employing a 
laser Doppler velocimeter, Frenklach et al. [55] measured the size of soot 
particles formed in the reactive flow. The particle size appeared to be of the 
order of a micron meaning that: 

a) agglomeration of the primary soot spheres takes place prior to the 
onset of cooling and 

b) the use of the Rayleigh limit of the Mie scattering theory for the 
determination of absolute soot yield is invalid and leads to an over- 
estimation of the conversion of hydrocarbon to soot. 

There are a number of current programs. Rawlins et al. [56,57] have been 
conducting an experimental program on the pyrolysis and oxidation of aromatic 
hydrocarbons behind incident shock waves. A variety of light extinction and 
emission measurements have been performed. Colket and Seery [58] have been 
conducting single-pulse shock-tube studies on the pyrolysis of toluene. Bauer 
and Zhang [59] have been performing single-pulse shock-tube studies on the 
pyrolysis of di- and tri-cyclic aromatics. 


b) Oxidation and Pyrolysis of Non-Aromatic Hydrocarbons in Shock Tubes . 


The chemistry of non-aromatic hydrocarbons is much better known than that 


of aromatics. Shock tube investigations have been recently reviewed by 
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Khandelwal and Skinner [60]. Since this review, a number of relevant papers 
have appeared. Frank and Just [61] reported decomposition of acetylene and 
diacetylene using atomic resonance absorption spectroscopy technique for 
hydrogen atoms. Bar-Nun and Dove [62] investigated acetylene pyrolysis and 
oxidation by water vapor. Koike and Morinaga [63] studied pyrolysis of 
acetylene and ethylene by UV absorption. Kern et al. [32] reported the results 
on pyrolysis of acetylene and 1,3-butadiene obtained behind reflected shock 
waves by time-of-flight mass spectrometry. Kiefer and co-workers [64-66] have 
studied the pyrolysis of ethylene, propane and propene by means of laser- 
schlieren technique. Colket [67] has been conducting experiments on the 
pyrolysis of 1,3-butadiene in single-pulse shock tube studies. 

The pyrolysis of acetylene has received the most attention. The 
experimental results obtained in a variety of kinetic studies were summarized by 
Tanzawa and Gardiner in a kinetic model [68,69]. The model reproduces the 
general sequence 

C 2 H 2 ♦ C 4 H 3 * C 4 H 2 * C 6 H 2 * C 8 H 2 ** * * 

suggested by Gay et al. [70] for acetylene pyrolysis and predicts the 
establishment of equilibrium among the smaller acetylenes and acetylenic 
radicals prior to soot formation. Modeling of acetylene oxidation has been 
recently reported by a number of researchers [71-73]. No reactions leading to 
soot, however, have been suggested. 

Graham et al. [43-45] were the first to monitor soot behind shock waves. 
The authors reported that much less soot was formed from acetylene than from 
aromatic hydrocarbons. Fussey and co-workers [74-76] employed a variety of 
optical techniques in studies of the pyrolysis of ethane, ethylene, acetylene 
and propylene. The authors determined that soot is formed more readily from 
acetylene than from ethylene and ethane. They determined the induction times 
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for soot appearance and postulated that these times are independent of 
experimental technique. Induction times for soot appearance in acetylene 
pyrolysis were also obtained by Tanzawa and Gardiner [69] and Yoshizawa et al. 
[77]. 

Frenklach et al. [78] have systematically investigated soot formation in 
pyrolysis of acetylene, allene and 1,3-butadiene. Soot was monitored by laser 
extinction at 632.8 nm and 3.39 pm behind reflected shock waves. The authors' 
major conclusions were the following: 

The induction times for soot appearance are dependent on the 
wavelength of extinction which invalidates the assumption of Fussey 
et al. [74] ; 

The formation of soot during acetylene pyrolysis is not only strongly 
dependent on the initial concentration of acetylene, the fact which 
has been reported by a number of researchers [42,69,74-77], but also 
exhibits the existence of a maximum in soot yield with temperature at 
higher concentrations of acetylene. The maximum occurs at 
significantly higher temperatures than in pyrolysis af aromatics; 
In acetylene pyrolysis, a decrease in total pressure shifts the soot 
bell to higher temperatures with a significant increase in the maximum 
soot yield; 

Soot is formed much faster and in much larger quantities from allene 
than from 1,3-butadiene; 

The most efficient "building blocks" for the formation of 
soot precursors in the pyrolysis of aliphatic hydrocarbons 
seem to be species with C-to-H ratios of approximately unity 
which have conjugated molecular configuration. Since aromatic 

r 

S 
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form provides the ultimate delocalization of rt-electron density 
and thus the ultimate stabilization and reactivity, the 
incipient soot formation from hydrocarbons should follow the 
route of consecutive production of the conjugated reactive 
structures. The difference in soot formation characteristics 
between various hydrocarbons is determined by the initiation 
process, i.e. by the reactions leading to these reactive 
structures . 

The conceptual model developed for soot formation is pyrolysis of non- 
aromatic hydrocarbons is discussed in Subsection 3b of DISCUSSION. 

2. Objective of the Present Work 

The main objective of the present work was to investigate soot formation 
characteristics in oxidation of aromatic and non-aromatic hydrocarbons. Toluene 
and acetylene were the chosen representatives of these two classes of 
hydrocarbons. These are the species for which soot formation dependencies were 
established in the previous pyrolysis studies [54,78]. The oxidation was 
conducted at the conditions of the previous pyrolysis experiments, in order to 
compare the results of both studies . 

It was recently reported by Wang et al. [46,47] and Rawlins et al. 
[56,57,80] that oxygen suppresses soot formation from toluene. The results 
obtained in the present work indicate that, depending upon experimental 
conditions, oxygen also promotes the production of soot. This also appeared to 
be true in the case of acetylene: depending on the experimental conditions, 

oxygen promotes or suppresses the formation of soot. 



17 


Some of the results obtained during this contract are already reported in 
the literature. The entire list of publications resulting from this contract is 
given in Appendix A. The professional personnel of LSU involved in this 
contract is given in Appendix B. 

APPARATUS AND PROCEDURES 


The experimental apparatus and procedures used in this study were similar 
to those described in our previous works on pyrolysis of aromatic hydrocarbons 
[57,78]. The experiments have been conducted behind reflected shock waves in a 
conventional stainless steel shock tube: 7.62 cm i.d., 3 m driver section and 
7.3 m driven section. The double diaphragm burst technique was employed, using 
either Mylar or aluminum foils as the diaphragm material. Both mechanical 
(Edwards ED-500) and diffusion (Edwards Speedivac E-04) pumps were used in the 
shock tube gas-handling and vacuum systems. The systems could be evacuated to 
less than 10 ^ torr. The driver gas was helium. The test gas mixtures were 
prepared manometrically in a stainless steel tank and allowed to mix for at 
least 24 hours prior to experimental runs. A Wallace and Tiernan Model FA-145 
manometer was used for precise pressure readings. The stated purities of the 
gases were: argon-99 .995%, helium-99 . 99% oxygen-99. 5% and acetylene-99 . 6%. The 
toluene (Reagent, Baker) and benzene (Spectranalyzed , Fisher) were purified by 
repeated freezing and evacuation. The shock tube was cleaned after every run. 

The state of the gas behind the reflected shock waves was calculated in a 
standard manner assuming full relaxation and no chemical reaction [79] and using 
the measured incident shock velocity extrapolated to the end wall of the shock 
tube. The observed shock wave attenuation was approximately 2%/m. Shock 
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velocities were measured using four piezoelectric pressure transducers to 
trigger the start and stop channels of an interval timer. All pressure 
transducers and optical windows were mounted flush with the surface of the shock 
tube to minimize flow distortion. 

The soot conversion was determined by measuring the attenuation of the beam 

from a 15 mw cw He-Ne Spectra-Physics laser which was operated in the visible 

(632.8 nm) region of the spectrum. The laser beam crossed the shock tube at a 

point 10 mm from the end plate. High quality sapphire windows were used. The 

attenuated laser light was monitored by an RCA 1P28 photomultiplier, whose 

output signal was amplified and displayed on a Nicolet Model 2090-3 digital 

oscilloscope. The design of the optical system was optimized so that emission 

was only a negligible component (less than 0.1%) of the extinction signal. This 

was achieved by using the laser beam at maximum power, a narrow-band 

interference filter centered at 632.8 nm and a number of optical stops. Careful 

alignment and adjustment of the optical system at 632.8 nm resulted in an 

excellent signal-to-noise ratio of the absorption signal. As in our previous 

studies [54,78], the induction time for soot appearance was defined by the point 

of maximum curvature in the extinction signal and the soot yields were 

calculated according to Graham's model [42,43] but leaving out the quantity E(m) 
2 2 

= - Im [(m -l)/(m +2)], where ra is the complex refractive index of soot 

particles. This arbitrary form of reporting the results was chosen to emphasize 
the ambiguity in the value of m and, what seems to be even more important, in 
the laser-extinction model itself [54,55]. 

In addition to the extinction measurements, the infrared emission was also 
monitored by employing two side-looking In-Sb detectors operated in a 
photoconductive mode. The detectors were located on the opposite side walls of 
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the shock tube, 65 mm away from the end plate. The fields of view were 
determined by vertical 1 x 4 mm slits. Two narrow-band interference filters, 
4.78±0.08 pm and 4.2510.05 pm, were used in an attempt to isolate the CO and CC^ 
emissions, respectively (the center of CO band is at 4.75 pm and the center of 
CO^ band is at 4.35 pm). 

The pressure behind reflected shock waves was monitored by a calibrated 
piezoelectric pressure transducer located on the upper wall of the shock tube, 
10 mm away from the end wall. Output signals from the pressure transducer, the 
1P28 photomultiplier and the two In-Sb detectors were amplified, displayed and 
recorded on Nicolet, Model 2090-3, digital oscilloscopes. The data were 
transferred from a Nicolet oscilloscope to an Apple II microcomputer via a 
GPIB/IEEE-488 interface. The data were reduced using PASCAL language programs 
in the Apple II before being transferred to the IBM main frame computer. 

RESULTS 

The initial conditions used in this experimental work are summarized in 
Table II. Mixtures A and C were used to test the reproducibility of the soot 
formation results, with data in literature. Mixtures B and D were used to study 
the effect of oxygen on soot formation from toluene. Mixtures E, F and G were 
used to complete the data on pyrolysis of toluene, which were used for empirical 
modeling of soot formation. Mixtures I, J, L, M and N were used to study soot 
formation in oxidation of acetylene. These experiments were designed so that 
the effects of temperature, pressure and oxygen concentration could be clearly 
observed. Mixtures 0 and P were used to study the pyrolysis and oxidation of 
benzene, respectively. 



20 


TABLE II. Experimental Conditions 


Mixture 

Composition 
(% vol. in argon) 
Fuel 0 2 

Temperature 

(K) 

Pressure 

(bar) 

C x 10 6 
(mol/cm 3 ) i 

[C] x 10" 17 
(carbon atoms/cm 3 ) 

A 

0.311, 

C 7 H 8 

- 

1525-2349 

1.98-3.00 

14.52-15.74 

1.90-2.06 

B 

0.311, 

C 7 H 8 

0.311 

1496-2391 

1.87-3.08 

15.03-15.49 

1.97-2.03 

C 

1.75, 

<*. 

- 

1684-2813 

0.32-0.66 

2.31- 3.13 

1.71-2.31 

D 

1.75, 

C 7 H 8 

1.75 

1668-2932 

0.33-0.68 

2.29- 3.15 

1.69-2.33 

E 

0.10, 

‘7*8 

- 

1647-2051 

2.07-2.70 

15.23-15.73 

0.64-0.67 

F 

0.5627 

’ C 7 H 8 

- 

1463-2480 

0.31-0.58 

2.54- 2.89 

0.60-0.69 

G 

1.0, 

C 7„8 

- 

1529-2932 

1.99-3.73 

14.83-16.10 

6.25-6.79 

H 

4.65 

C 2 H 2 

- 

1687-3123 

1.25-2.33 

8.70- 9.22 

4.87-5.16 

I 

4.65 , 

C 2 H 2 

4.65 

1532-3490 

1.10-3.14 

7.91-11.05 

4.43-6.19 

J 

4.65, 

C 2 H 2 

1.5 

1429-2259 

1.08-1.70 

8.74- 9.34 

4.89-5.23 

K 

20.00 

C 2 H 2 

- 

1748-2802 

0.28-0.49 

1.95- 2.25 

4.69-5.41 

L 

20.00, 

C 2 H 2 

6 . 45 

1379-2623 

0.22-0.37 

1.73- 2.48 

4.16-5.96 

M 

20.00, 

c 2 h 2 

4.30 

1431-2758 

0.20-0.40 

1.67- 2.42 

4.01-5.83 

N 

20.00, 

C 2 H 2 

2.15 

1578-2997 

0.26-0.50 

1.70- 2.36 

4.09-5.68 

0 

0.311, 

‘6*6 

- 

1561-2273 

1.94-2.87 

14.90-15.59 

1.67-1.75 

P 

0.311, 

C 6*6 

0.311 

1516-2243 

1.98-2.83 

15.16-15.68 

1.70-1.76 


These ’experiments were performed with the partial support of the Office of 
Fossil Energy, U.S. Department of Energy, under the auspices of Grant Number 


DE-FG22-80PC30247. 
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Figures 1 and 2 present typical experimental records and Figure 3 explains 
our definitions of the induction time for soot appearance, t soot > and the rate 
of soot production, R soot - There were some instances when the induction time 
and/or rate of soot formation could not be measured. They were: 1) when the 
shape of the laser extinction trace was such that the inflection point on the 
trace could not be clearly defined, which was particularly true for acetylene 
mixtures at low pressures; 2) at very low temperatures when ' t soot was longer 
than the observation time; 3) at very high temperatures when the induction time 
was very small or close to zero. It might still be possible to measure the rate 
of soot formation under condition 3). 

Comparing the corresponding infrared emissions obtained at similar 

conditions in pyrolysis and oxidation of toluene (Figs. 1 and 2), we noted a 
slight decrease in the emission when oxygen was added. Since we would expect no 
CO or CO^ to be present in pyrolysis, the observed decrease in the emission 
signals indicates that these signals are primarily due to the "black body 
radiation" of soot particles and the contribution from CO and CO^ molecules is 
probably negligible. Therefore, we must conclude that it is impossible to 
monitor CO and CO^ under the conditions employed. Rawlins and co-workers [80] 
came to a similar conclusion in their work. Therefore, no further attempts to 
monitor CO and CO^ were undertaken in the rest of the study. 

The measured soot yields, induction times for soot appearance and rates of 
soot formation are reported in Figures 4 through 25. Figures 4-7 show the 
effect of oxygen on soot formation from toluene at high pressures. As can be 
seen from these figures, oxygen suppresses the formation of soot; this effect is 
qualitatively somewhat similar to that resulting from a reduction in the initial 
concentration of toluene (see Fig. 26). At lower pressures (Figs. 8-11), the 
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FIGURE 7. Rates of soot formation versus inverse temperature in 
toluene mixtures at high pressures. 
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The influence of oyxgen on soot yield in toluene 
mixtures at low pressures. 
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FIGURE 13. The influence of oxygen on soot yields in acetylene 
mixtures at high pressures. 


Logio ( T S oot ) » (s 0 c) 


ABSORPTION AT 632.8 nm 
CC3-N-5 . 00 *10 17 carbon atoms/cm 3 

^ 4.65 XC 2 H 2 + Ar 

P=1 ; 25 t o 2 . 33 bar 

1- 4.65 7.C 2 H 2 + 1 . 50 7.a 2 +Rr 

P=1 . 08 to 1.70 bar 


lO 4 /T, 


FIGURE 14. Induction times for soot appearance versus inverse 
temperature in acetylene mixtures at high pressures 
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FIGURE 15. Rates of soot formation versus inverse temperature in 
acetylene mixtures at high pressures. 
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FIGURE 16. Soot yields versus temperature at different reaction 

times in 20.0% acetylene - 4.30% oxygen-argon mixtures. 
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FIGURE 17. Soot yields versus temperature at different reaction 

times in 20.0% acetylene - 2.15% oxygen-argon mixtures 
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FIGURE 19. Soot yields versus temperature at different reaction 
times in 0.311% benzene-argon mixtures. 
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FIGURE 20. Soot yields versus temperature at different reaction 

times in 0.311% benzene - 0.311% oxygen-argon mixtures. 
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FIGURE 21. Comparsion of soot yields obtained in toluene and 
benzene pyrolysis. 
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FIGURE 22. Comparison of soot yields obtained in toluene and benzene 
oxidation. 
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FIGURE 23. The influence of oxygen on soot yield in benzene mixtures. 
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FIGURE 24. Comparison of induction times for soot appearance in 
toluene and benzene mixtures. 
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FIGURE 25. Comparison of rates of soot appearance in toluene 
and benzene mixtures. 
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FIGURE 26. Soot yields versus temperature in toluene-argon mixtures 
at a reaction time' of 1 ms. 
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addition of oxygen not only reduces the maximum soot yield but also shifts the 
maximum to a lower temperature. In other words, at low pressures and low 
temperatures oxygen actually promotes soot formation. 

Figures 12 through 18 present the results obtained with acetylene mixtures. 
In larger amounts, oxygen completely suppresses the production of soot (Mixtures 
I and L) . The addition of oxygen in smaller quantities shifts the soot yield 
"bell" to lower temperatures. This shift is much more pronounced than that in 
the toluene case. Furthermore, the shift in the case of acetylene is observed 
throughout the pressure range tested while the shift was observed in toluene 
mixtures only at lower pressures. It is interesting to note, however, that the 
sensitivity of the soot yields to oxygen concentration in acetylene mixtures is 
much higher at lower pressures than at higher ones. 

Figures 19- 25 present soot yields observed in pyrolysis and oxidation of 
benzene. As can be seen in these Figures, there- is no significant difference 
between the results obtained in benzene and toluene mixtures . 

Figure 26 presents the results obtained in pyrolysis of toluene. The 
conditions used in these mixtures allow one to determine the dependence of soot 
yield on the initial concentration of toluene. The reasons for this dependance 
can be qualitatively understood by referring to a recently published conceptual 
model [54] . 

DISCUSSION 


The most striking experimental observation is the shift of the soot-yield 
bell to lower temperatures when oxygen is added. However, prior to discussion 
of this result, one would like to estimate the temperature change during the 
reaction. The overall combustion mechanism is composed of many elementary 
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reactions, each either endothermic or exothermic; therefore, it is most probable 
that the temperature behind the reflected shock wave will change with time. 
Since we cannot measure the temperature in situ , it was decided to estimate the 
temperature change, which may occur in combustion of toluene and acetylene under 
our experimental conditions, via computer modeling. 

1. Estimation of Temperature Change During Oxidation of Toluene . 

To simulate the oxidation of toluene, we chose the mechanism of Jachimowski 
[81-83]. The mechanism involves 24 chemical species and 52 elementary 
reactions. The reactions and their forward rate constants are given in Table 
III. The reverse rate constants were calculated using the equilibrium 
constants, the thermo chemical data for which were taken from the updated NASA 
polynomials of Burcat [84], 

It is not suggested here that Jachimowski ' s mechanism represents the 
reaction process exactly; for example, it does not account for soot formation. 
However, no other mechanism is available which is clearly superior for our 
purposes. Although better estimates of some of the rate constants are now 
available, Jachimowski presented a unified mechanism which had rate constants 
adjusted to reproduce observed experimental data obtained in shock tubes. We 
believe, therefore, that the temperature changes predicted by the Jachimowski 
mechanism can be expected to match actual values reasonably well, allowing the 
estimation of possible variations in the temperature behind reflected shock 
waves at our conditions . 

The Jachimowski mechanism was incorporated into a constant-density computer 
model which used the latest version (1982) of LSODE [85]. The gas mixtures used 
in the simulation were those of Mixture B, i.e. 0.311% toluene - 0.311% oxygen 
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-argon, and Mixture D, i.e. 1-75% toluene - 1.75% oxygen-argon. Species 
concentrations and temperature profiles were calculated to a reaction time of 
3000 microseconds for five initial temperatures: 1600K, 1800K, 2000K, 2200K, 
and 2400K. These temperatures were chosen as typical of the operating 
temperature range of interest. The temperatures at various reaction times 
calculated for all five cases are summarized in Tables IV and V. 

Inspection of the simulation results indicates that the temperature change 
as the result of chemical reaction in the first three milliseconds for Mixture B 
is not significant, being on the order of magnitude of the uncertainty in the 
initial post-shock temperature, i.e. ~25K. As expected, somewhat larger 
temperatures changes are seen for Mixture D. However, the kinetic simulations 
showed that the maximum change in temperature was still moderate, being about 
80°K. The typical temperature versus reaction time profile appears to be a 
shallow ascending s-curve. The initial temperature decrease is caused by the 
following endothermic reactions: the attack of on C^Hg to form the free 
radicals C^H^ and HO^ (reaction 1) and the fragmentation of toluene (reaction 
2). As the reaction time increases, exothermic reactions involving oxygen 
become important. These reactions form relatively stable products such as ^0, 
CO, and CC^. However, endothermic decompositions are occurring simultaneously, 
moderating the rise in temperature. At longer times, when most of the oxygen is 
bound in relatively stable products, the continuing pyrolytic fragmentation of 
the remaining toluene causes the final temperature decrease. 

We also performed several adiabatic equilibrium calculations using a 
computer program developed at the NASA Lewis Research Center [86]. The results 
of these computations are given in the last column of Tables IV and V. The 
adiabatic equilibrium temperature is expected to provide an upper limit for 
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TABLE IV 

Temperatures Computed for the 0.311% Toluene 

- 0.311% Oxygen- Argon Mixture at a carbon atom 

17 3 

concentration of 2.0 x 10 atoms/cm . 


Initial 

Temperature 

(K) 

Temperature Computed by Kinetic Model 
at Reaction Times, in ps 

Equilibrium 

Temperature 

(K) 

100 

500 

1000 

3000 

1600 

1590 

1588 

1596 

1617 

1689 

1800 

1779 

1809 

1814 

1798 

1892 

2000 

1982 

1989 

1983 

1978 

2094 

2200 

2192 

2197 

2197 

2197 

2292 

2400 

2384 

2385 

2385 

2387 

2483 


TABLE V 


Temperatures Computed for the 1.75% Toluene - 

1.75% Oxygen-Argon Mixture at a carbon atom 

17 3 

concentration of 2.0 x 10 atoms/cm . 


Initial 

Temperature 

(K) 

Temperature Computed by Kinetic Model 
at Reaction Times, in |Js 

Equilibrium 

Temperature 

(K) 

100 

500 

1000 

3000 

1600 

1565 

1554 

1576 

1675 

1996 

1800 

1717 

1819 

1845 

1783 

2195 

2000 

1913 

1964 

1934 

1904 

2375 

2200 

2157 

2155 

2155 

2154 

2532 

2400 

2347 

2373 

2382 

2389 

2668 
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temperature change since it assumes complete conversion to the stable products, 
whether or not production of these species is kinetically possible in the 
available reaction time. These equilibrium calculations, however, do not 
provide a precise limit, particularly in fuel-rich combustion, because many 
complex hydrocarbon species are not in the computer program data base. The 
difference between the initial post-shock temperature and the calculated 
adiabatic equilibrium temperature for Mixture B was moderate, ranging from 83°K 
for an initial temperature of 2400K to 94°K for an initial temperature of 2400K. 
The difference between the initial temperature and the adiabatic equilibrium 
temperature was, as expected, much larger for Mixture D. 

In summary, the conclusion drawn from our computer modeling of toluene 
oxidation is that for our experimental conditions and available reaction time, 
the temperature change as a result of chemical reaction is not significant 
compared to the temperature shift of the soot bell observed experimentally. 
This conclusion is supported for most of the mixtures by the results of 
adiabatic equilibrium calculations. 

2. Estimation of Temperature Change During Oxidation of Acetylene 

To estimate the temperature change during oxidation of acetylene at the 
conditions used in this study, we chose the mechanism proposed by Gardiner and 
co-workers [71,87,88], the most recent mechanism known to us. 

The mechanism is based in part on work done by other researchers at various 
experimental conditions, but it was "tuned" to match their own recent work with 
fuel-lean mixtures (cj) = 0.25 to 1.0). This mechanism therefore does not exactly 
represent the reaction process for our fuel-rich conditions (<[> = 7.75). However, 
we believe that the temperature changes predicted by Gardiner's mechanism can be 
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expected to predict the actual values within at least an order of magnitude, 
allowing the rough estimation of possible temperature variations behind 
reflected shock waves at our conditions. 

The computer program is a constant-density model which uses the latest 

version (1982) of LSODE. The computer simulations were performed for the 

conditions of Mixture J (see Table II), i.e. 4.65% acetylene - 1.5% oxygen - 

argon, and Mixture M, i.e. 20.0% acetylene - 4.3% oxygen-argon. Both mixtures 

had a carbon atom concentration behind the reflected shock of 5.0 x lO 1 ^ 
3 

atoms/cm . Species concentrations and temperature profiles were calculated for 
five initial temperatures: 1500K, 1600K, 1700K, 1800K, and 1900K for Mixture J 

and 1600K, 1800K, 2000K, 2200K, and 2400K JJor—Mixture M. These temperatures 
were chosen as typical of the operating temperature ranges of interest (see 
Figs. 12 and 16). The temperatures at various reaction times calculated for all 
cases are summarized in Tables VI and VII. 

Inspection of the computational results indicates that the temperature 
change as a result of chemical reaction, for our gas mixture, experimental 
conditions, and reaction time, is not significant, being less than 50K. The 
small temperature rises seen are due to the quick exothermic formation of such 
relatively stable species as CO, anc * ^O. ® nce raost of the available 

oxygen is bound in these products, the remaining acetylene reacts very slowly, 
causing little temperature change. 

Somewhat larger temperature changes were predicted by adiabatic equilibrium 
calculations, especially at low initial temperatures. This was due to the fact 
that the NASA computer program [86] used for the equilibrium calculations 



TABLE VI. 


Temperatures Computed for the 4.65% Acetylene - 

1.5% Oxygen-Argon Mixture at a carbon atom concentration of 

17 3 

5.0 x 10 atoms/cm . 


Initial Temperature Computed by Kinetic Model Equilibrium 

Temperature at Reaction Times in (Js Temperature 


(K) 

50 

250 

500 

2500 

(K) 

1500 

1500 

1502 

1501 

1501 

2296 

1600 

1600 

1602 

1602 

1602 

2363 

1700 

1700 

1703 

1702 

1702 

2422 

1800 

1803 

1804 

1803 

1804 

2472 

1900 

1905 

1906 

1905 

1906 

2514 
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TABLE VII. 

Temperatures Computed for the 20% Acetylene - 

4.3% Oxygen-Argon Mixture at a carbon atom 

17 3 

concentration of 5.0 x 10 atoms/cm . 


Initial Temperature Computed by Kinetic Model Equilibrium 

Temperature at Reaction Times in (Js Temperature 


(K) 

50 

250 

500 

2500 

(K) 

1600 

1600 

1600 

1601 

1599 

3107 

1800 

1800 

1806 

1806 

1804 

3167 

2000 

2008 

2012 

2012 

2012 

3225 

2200 

2215 

2221 

2224 

2228 

3282 

2400 

2424 

2435 

2441 

2443 

3339 
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includes soot-like species such as C^qH-jq as possible reaction products. The 
equilibrium calculations indicated the production of significant amounts of 
these large soot-like species. However, at our reaction times, the system is 
still in the kinetic regime for soot and is far removed from complete 
equilibrium. 

In summary, the conclusion drawn from our computer modeling of acetylene 
oxidation is that for the experimental conditions tested and the available 
reaction times, temperature change as a result of chemical reaction is not 
significant compared to the experimental temperature shift of the soot bell 
under oxidative pyrolysis as compared to pure thermal pyrolysis. 

3. Soot Formation in Hydrocarbon Oxidation 

a) Toluene Mixtures 

Comparing the experimental results on soot formation obtained in pyrolysis 
and oxidation of toluene, we conclude that while the addition of moderate 
amounts of oxygen can have an effect on soot levels and temperature regions 
where soot is formed, it seems to have little effect on the general character of 
soot formation observed in pyrolysis. The time development of the soot "bell" 
and its shift to higher temperatures when pressure is lowered are similar in 
both cases. Actually, addition of oxygen to toluene seems to be qualitatively 
equivalent to reduction of the initial concentration of toluene (cf. Figs. 5 
& 26). 

As in the pyrolysis case, pressure has little effect on the maximum soot 
yield values in oxidation. However, the observed shift of the soot yield 
maximum is smaller when oxygen is present (cf. Figs. 5 & 9). In other words, at 
relatively low pressures the addition of oxygen results in the shift of the soot 
yield maximum to lower temperatures while at relatively high pressures the shift 
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is hardly observable. This observation can be also interpreted as that at low 
pressures oxygen suppresses soot production at higher temperatures whereas at 
lower temperatures oxygen promotes the formation of soot (Fig. 9). 

The observed phenomenon indicates that addition of oxygen causes reactions 
to occur that compete with pressure-dependent processes. In our recent analysis 
of soot formation during pyrolysis of aromatic hydrocarbons [54] , we suggested 
that pressure-dependent fragmentation of an aromatic ring initiates the 
production of soot. Hence, oxidation reactions must compete with ring 
fragmentation. At low pressures, the promoting effect of oxygen on soot 
formation below 2100K and the suppression of soot above this temperature 
indicate that oxidizing agents (such as 0, OH) attack the aromatic ring 
producing, besides oxidation products CO and CO2 and small inefficient 
fragments, active intermediates for the soot-forming route. At low temperatures 
when thermal ring fragmentation is slow, the formation of these intermediates by 
oxidative reactions enhances soot production. At high temperatures, when the 
thermal decomposition is predominant, additional removal of aromatic rings by 
oxidative agents inhibits formation of soot. 

b) Acetylene Mixtures 

In the case of acetylene (Figs. 13 and 18), the effect of oxygen on soot 
yield is more pronounced and extends over the entire pressure range studied. It 
is interesting to note that the sensitivity of the acetylene system is much 
larger than that of the toluene system; that is, a smaller oxygen- to-hydrocarbon 
ratio is required to produce a similar effect. The sensitivity to oxygen in the 
acetylene case is particularly high at lower pressures (Fig. 18). 

Soot formation in pyrolysis of acetylene was explained recently [78] in 
terms of the following conceptual model: 
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C 0 H 0 -> X 

2 2 ~ 


where: 




A 



X 



S 


M = a third body 

X = {X^, X 2 >...} is the collection of nonaromatic intermediates; 

A = {A^, A^,...} is the collection of aromatic species; 

s = {S S is the collection of species absorbing the light at 

a specific wavelength. 

The products of acetylene pyrolysis, X, interact among themselves, eventually 
forming the aromatic species A. The reactions of nascent aromatic species A 
comprise two parallel routes: the pressure-dependent high-activation-energy 

fragmentation of an aromatic ring and a low-activation-energy radical-molecule 
interaction of an aromatic ring with aliphatic fragments leading to soot. 

At low temperatures, the rate limiting step in acetylene pyrolysis is the 
production of X 

C 0 H_ X. 

Addition of relatively small amounts of oxygen results in formation of reactive 
radicals (e.g. 0 and OH), which promote the pyrolysis reactions. Under these 
conditions, the rate of formation of aromatic intermediates A is enhanced, 
which, in turn, increases the rate of soot production. 

At high temperatures, a partial equilibrium is quickly established and the 
fragmentation of aromatic rings 
A -> X 

<v /v 

becomes a dominant factor in soot formation process. As was postulated earlier, 
reactions of the oxidizing agents with A results in additional removal of 
aromatic rings which inhibits the formation of soot. 



63 


Thus, the overall effect of oxygen on soot formation in acetylene pyrolysis 
should be the shift of the soot yields to lower temperature which was observed 
experimentally. The higher sensitivity of soot production to oxygen for 
acetylene compared to that for toluene can also be explained now. At low 
temperatures, oxygen very efficiently promotes the formation of aromatics from 
acetylene compared to a relatively marginal contribution to the formation of X 
in toluene case. At high temperatures, oxygen promotes the ring fragmentation, 
which effectively retards the production of soot for both cases. 

c) Benzene Mixtures 

Soot formation in pyrolysis and oxidation of benzene is very similar, both 
qualitatively and quantitatively, to that of toluene. This fact provides 
further support for our conceptual model for soot formation in pyrolysis and 
oxidation of aromatic hydrocarbons. 

d) General Conclusions 

The main conclusion to be drawn from the results of this work is that the 
soot formation mechanism is probably the same for both pyrolysis and oxidation 
of hydrocarbons. That is, the addition of oxygen does not alter the soot route 
but rather promotes or inhibits this route by means of competitive reactions. 

The above conclusion actually implies that radicals are more important than 
ions in the soot formation mechanism since one would expect many more ions to be 
present during oxidation than during pyrolysis. This is particularly obvious in 
the case of aromatics (toluene and benzene). If ions were the crucial 
intermediates for soot precursors, one would expect a very dramatic positive 
effect on soot yields when oxygen is added. On the contrary, the observed 
effect is relatively small and negative. 
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4. Empirical Modeling of Soot Formation 

Quantitative prediction of the amount of soot formed in practical 
combustion systems has become one of the foremost concerns in combustion science 
and technology. Hitherto, no attempt has been made to develop a quantitative 
model, physical or empirical, that will relate the amount of soot formed with 
experimental conditions such as temperature, pressure, fuel concentration, et 
cetera. Previously, we have introduced a conceptual model for soot formation in 
pyrolysis of hydrocarbons which provides a starting point for empirical 
modeling. 

An approach to empirical modeling of soot formation is presented below. 
The objective of the modeling is to predict soot yield for a given reaction time 
at various temperatures, pressures and initial reactant concentrations for 
various fuels. The approach is introduced by first considering the simplest 
case, which is soot formation in shock-tube pyrolysis of aromatic hydrocarbons. 


a) Statement of the Problem 

The empirical model is postulated as 


A + X Q 


A + X, 


A + X.- 
1 


- X i + 1 


fragmentation 


polymerization 


Cl) 


which is the simplest, in a mathematical sense, kinetic scheme that has all the 
features of our conceptual model for soot formation in pyrolysis of aromatic 
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hydrocarbons [54].- Species A denotes aromatic hydrocarbon and/or their stable 
radicals and species X^, X K+ ^, ... of this system constitute what we will call 

soot, i.e. 

s = iv X K +1 . •••!■ (2) 

The soot yield is then defined as 
00 

Y= I (i+l){[X ]/[A] } (3) 

i=K 

where: Y is the soot yield; 

[X.] is the concentration of species X^; 

[A]g is the initial concentration of A. 

"Yield" is thus the fraction of species A originally present which has 

become soot. The (i+1) factor accounts for the number of A molecules 

"contained" in each X. molecule. 

l 

Let us derive an expression for Y as function of experimental parameters. 


b) Derivation 

The differential equations for the reaction system (1) take the form of 


d[A] _ 
dt 

d[X Q ] 

dt 

d[X 1 ] 
“dt = 


k f [A] - ,| 0 k p [A][X.] 


k f [A] - k p [A][X Q ] 


k P t A H x 0 ] - y A nv 


d[ x i+ l] 

-dT^ = kp[A][X.] - k p [A][X. +1 ] 


( 4 ) 
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with initial conditions 

^ t =0 - £ a ] 0 

[X.] t=0 = 0 , i = 0 , 1 , 2 ,... . 

Introducing the "dimensionless concentrations" 
a = [A]/[A] Q and x ± = [X ± ]/ [A3 q , 

dividing both sides of each differential equation (4) by k^a, and substituting, 
according to [89], 

k^adt = dx (5) 

we obtain 


da 

dt 


00 

1 - .2 n x./r 
J-0 J 



(6a) 


(6b) 




67 


P(i+1 


T/r 

, T/r) = I e U u 1 du 

* / 0 


is the incomplete Gamma function [90]. It can be noted that 


00 


.1 x. = T . 

1=0 l 


Indeed, 


oo 


00 


.1 x. = .1 r P(i+1, t/r) = r .I n — 

i=0 i i=0 i=0 l! 

T/r T/r 


T/r 

f-u . 

I e u du 


)du 


as /* -u i /* “ -u i 

■ ■ .i. J Hr V- 
"7 

= r I du = T, 


since by definition of the Poisson distribution [91] 

00 -u i 

1 £— g- = 1 

i=0 i! X * 

Substituting (10) into (6a), we obtain 
dT ~ “ 1 X/r ’ 

the solution of which, taking into account the initial conditions 
a = 1 - T - T 2 /2r 

Substituting (14) into (5), we obtain 
t T 

j ~ k f dt = J" dt/ (1 - T - T 2 /2r) 

*0 0 
or 

T/A 0 -r 


1 “'2 

k t = - In ( — jz ) 

f q X/X -r 


(9) 

( 10 ) 


(ID 

( 12 ) 

(13) 

(6c), is 

(14) 



68 


where: 

\ = 'i + q; 

\ 2 = -l - q; 
q = V 1 + 2/r; 
or 

qk f t 

T = rd-e ) 

1 1 qk f t 

\ 2 " \ e 

The last expression relates the "transformed time" T, the physical meaning of 
which is disclosed through equality (10), with the real time t. Hence, 

substitution of (7) and (15) into (8) and (14) determines the kinetic behavior 
of reaction system (1). 

The expression for soot yield (3) can be developed in the following manner. 
Let us rewrite expression (3) as 

oo oo k-1 

Y = I (i+l)x = I (i+l)x. - I (i+l)x. (16) 

i-K 1 i=0 1 i=0 x 

where the first sum on the right-hand side of this expression is easily 
determined: 


00 


00 


00 


I (i+l)x. = I ix. + I x.. 
i=0 1 i=0 1 i=0 x 


(17) 


However, I x., according to (10), equals I and, as in derivation (11), we 
i=0 X 


obtain 


T/r 

i-0 ix i = ill i,r P(i+1 - x /r) = * ih i ’ h „/ e ' V 


du 


X/r X/r 

Oo / - u i / Oo - u l 

- r A i frfit = 


du 
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T/r X/r 

/ “ -u j / 

= r / u ( I e 77 ) du = r / udu = x 2 /2r 


( 18 ) 


where : 


00 -u j 

. _ . , 3 ^ e u J 

J - l ' 1 and j— o — 


= 1 by the definition of Poisson's distribution 
[91]. 


That is 


00 


I (i+l)x. = t + t 2 /2r 
i=0 1 


(19) 


In order to determine the second sum on the right-hand side of expression 
(16), let us recall one of the properties of the incomplete Gamma function [90], 
i . e . 


-t/r, . .i 

P(i+l,T/r) = P(i,T/r) - i ; ~- 


( 20 ) 


or 


P(i+1, x/r) = 


, | e- t/c (T/r) J 

j=0 j! 


( 21 ) 


Therefore , 


x . 

i 


1 e' X/r ft/rl j 

' t 1 - jio ; j[ x/r) 1 


( 22 ) 


and then 


K-i K-l i -x/r ( 

ilo f 1 - jio — ? J) > 


K-l K-l i -x/r, . .j 

r [ iioU^D - l Ci*» -l g ] 


K-l 


= r [ 


k(k+ 1) "■ * (i+l+K) (K-i) e T ^ r (x/r) 1 


i-0 


1 ! 


(23) 


Thus, substituting (19) and (23) into (16), we obtain 
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Y = r [t/r + l 2 /2r 2 


K(K+1) K : 1 (i+l+K)(K-i) e' T/r (l/r) i , 

2 i-0 2 i ! J 


c) Model Fit 

Expression (24) together with relationships (7) and (15) constitute an 
empirical model for soot formation, which contains 3 unknowns: K, k^, and k^. 

Assuming Arrhenius form for the rate constants 
k f = A exp (-6 f /T)p a 
and k p = B exp (-6 p /T)p^[A]^ 
where : 


T is the initial reaction temperature in K; 

3 

p is the total density in mol/cm ; 

3 

and [A] Q is the initial concentration of toluene in mol/cm , 
there are 8 parameters to determine, namely A,B,0j,8p,a,p,y and K. Although the 
quality of the fit was slightly better for K=3, for physical reasons [59] we 
assumed k=6. 

A general approach for fitting the model would be to optimize 


2 (Y , . - Y . . J , 

^ t calc,i,t expt,i,t ’ 


(25) 


where Y . . and Y , . are the experimentally observed and the calculated 
soot yields for i-th experiment at the observation time t. The summation is 
taken over the entire set of experiments and the chosen number (25 in this work) 
of time intervals. However, the experimental determination of soot yields 
strongly depends on the knowledge of the complex refractive index m[54]. Since 
the value of m is not well-established, it was desirable that the modeling 
results would not depend on m. This was achieved as follows. Instead of 
optimizing the objective function (25), we minimized 
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where is the experimental soot yield at specified, reference conditions 

and Y* , is the soot yield calculated by the model at these conditions. The 
calc 3 3 

reference conditions in our modeling are the 0.311% C^H 0 -Ar mixture at T = 1977 

7 o 


-5 3 

K, p = 1.54 x 10 mol/cm and time = 1.0 ms. The conditions corresponded to an 
actual experimental point taken at the approximate center of the ranges of the 
experimental variables of interest. 

Actually, since the present modeling is empirical in nature and one should 
not expect that Eq. (24) predicts the true values of soot yields, the objective 


function for optimization (26) should be rewritten as 


(27) 


The form of expression (27) implies that (Y/Y*) ^ should be considered as a 

single entity and no physical meaning should be attached to Y ^ or Y* ^ 
alone. That is, the empirical model predicts the relative values and in order 
to obtain the absolute value of soot yield at given conditions one must multiply 
the corresponding (Y/Y*) ca ^ c by Y* x p t . 

Optimizing (27) for all 5 toluene mixtures with no oxygen present, A,C,E,F 
and G, the following results were obtained: 


A 

— 

6.67 x 

10 

9 f 

= 

2.38 x 

10 

a 

= 

0.858 


B 

= 

2.68 x 

10 

6 

P 

= 

6.47 x 

10 

P 

= 

0.139 


Y 

= 

-0.413 


K 

— 

6 
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Figures 27-29 compare the experimental and computed values of Y/Y*. Figure 27 

shows the time development, Fig. 28 - pressure dependence, and Fig. 29 - 

concentration dependence of soot yields. Although the quality of fit is 

generally good, it can be improved if Mixture G is not taken into account. 

Thus, by optimizing (27) for only 4 mixtures, A,C,E and F, the following results 

were obtained: 

A = 1.94 x 10 14 

e f = 2.68 x 10 4 

a = 1.27 

B = 6.05 x 10 8 

0 = -4.38 x 10 3 

P 

p = -0.361 

Y = 0.0114 

K = 6 

Figures 30 - 32 compare the experimental and new computed values of Y/Y*. 
As can be seen in these figures the quality of fit is significantly improved. It 
is not clear, however, whether this improvement is due to some experimental 
problems with Mixture G or simply the result of the empirical nature of the 
modeling. 

Although it is very tempting to assign physical meaning to the rate 
constants above, one should be very cautious about doing so. These constants 
are parameters of an empirical model. 
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FIGURE 27. Comparison of the model prediction with experimental 

observations. (First parameter set). Time development. 
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FIGURE 29. Comparison of the model prediction with experimental 
observations. (First parameter set). Concentration 
depend ence. 
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FIGURE 30. Comparison of the model prediction with experimental 
observations. (Second parameter set). Time 
development . 
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FIGURE 31. Comparison of the model prediction with experimental 
observations. (Second parameter set) . The effect of 
pressure. 
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FIGURE 32. Comparison of the model prediction with experimental 
observations. (Second parameter set). Concentration 
dependence. 



SUMMARY OF RESULTS 


1. A conceptual model for soot formation in pyrolysis of non-aromatic 
hydrocarbons was developed. The most efficient "building blocks" for the 
formation of soot precursors in the pyrolysis of aliphatic hydrocarbons seem to 
be species with C-to-H ratios of approximately unity which have conjugated 
molecular configuration. Since aromatic form provides the ultimate 
delocalization of n-electron density and thus the ultimate stabilization and 
reactivity, the incipient soot formation from hydrocarbons should follow the 
route of consecutive production of the conjugated reactive structures. The 
difference in soot formation characteristics between various hydrocarbons is 
determined by the initiation process, i.e. by the reactions leading to these 
reactive structures. 

2. Soot formation in toluene-, benzene-, and acetylene-oxygen-argon mixtures 
was investigated to study soot formation in a combustion environment. It was 
observed that high concentrations of oxygen completely suppress soot formation. 
The addition of oxygen at relatively low concentrations uniformly suppresses 
soot formation at higher pressures, while at relatively lower pressures it 
suppresses soot formation at higher temperatures while promoting soot production 
at lower temperatures. The observed behavior may indicate that oxidation 
reactions compete with ring fragmentation. The main conclusion to be drawn from 
the results of this work is that the mechanism of soot formation in shock tubes 
is probably the same for both pyrolysis and oxidation of hydrocarbons. That is, 
the addition of oxygen does not alter the soot route but rather promotes or 
inhibits this route by means of competitive reactions. The above conclusion 
actually implies a radical and not an ionic mechanism of incipient soot 
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formation in shock tubes (although nascent ions may play an important role in 
coagulation). Indeed, if ions were the crucial intermediates for soot 
precursors, one would expect a very strong increase in soot yields when oxygen 
is added. On the contrary, the observed effect is relatively small and 
negative. 

3. An approach to empirical modeling of soot formation is suggested. An 
empirical model for soot formation in shock-tube pyrolysis of aromatic 
hydrocarbons is developed. The model predicts well the times, concentration and 
pressure dependecies of experimental soot yields. 
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FUTURE RESEARCH DIRECTIONS 


The ultimate goal of our program is the development of a detailed kinetic 
model for soot formation in combustion of hydrocarbons. To achieve this goal, 
we will concentrate on the following objectives: 

1 . Systematic experimental investigation of sooting characteristics at care- 
fully designed conditions . The conditions will be chosen to help in the 
identification of the essential steps in the mechanism of soot formation 
and also to provide sufficient information for empirical (part 2) and 
kinetic (part 3) modeling. This part of the program will be directed 
towards the establishment of the temperature, pressure, initial hydrocarbon 
concentration and molecular structure effects on soot yield for carefully 
selected hydrocarbons, heterocyclic compounds (for example, pyridine and 
chlorobenzene) and their mixtures. Both pyrolysis and oxidation will be 
performed. The experiments will be conducted behind reflected shock waves; 
soot will be monitored by the attenuation of a laser beam. Additional 
optical diagnostics will be employed if it will conclusively benefit the 
stated objectives. 

2. Empirical modeling of soot formation . In this part of the program empirical 
relationships between soot yield and temperature, pressure and composition 
of the mixture will be developed based on the results obtained in part 1. 
The relationships will be tested for flames. The established correlations 
may be used, on one hand, for practical purposes such as design and 
optimization of combustors, and, on the other hand, will be imposed as 
constraints for detailed kinetic modeling in part 3 of the program. 
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3. Detailed kinetic modeling of soot formation . Established mechanisms for 
pyrolysis and oxidation of hydrocarbons will provide a basis for the 
modeling. The reactions leading to the formation of soot precursors and 
their subsequent growth will be guessed based on analysis of the results 
obtained in part 1 and those of other researchers. The missing 

thermochemical data will be estimated. The kinetic model will be subjected 
to constraints established in parts 1 and 2 of the program and taken from 
other works. Sensitivity analysis will be employed to identify the main 


reaction route. The model will be extended to flames. 
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